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Abstract 

The aim of this paper is to derive Agmon's type exponential estimates 
for solutions of elliptic systems of partial differential equations on R n . We 
show that these estimates are related with the essential spectra of a family 
of associated differential operators which depend on the original operator, 
and with exponential weights which describe the decrease of solutions at 
infinity. The essential spectra of the involved operators are described by 
means of their limit operators. 

The obtained results are applied to study the problem of exponential 
decay of eigenfunctions of matrix Schrodinger, Moisil-Theodorescu, and 
Dirac operators. 



1 Introduction 

The main aim of the paper is to obtain Agmon's type exponential estimates 
for the decaying behavior of solutions of systems of elliptic partial differential 
equations with variable coefficients. Exponential decay estimates of this type 
are intensively studied in the literature. We only mention Agmon's pioneering 
papers [TJ [2] where estimates of eigenfunctions of second order elliptic operators 
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were obtained in terms of a special metric, now called the Agmon metric. Ex- 
ponential estimates for solutions of pseudodifferential equations on R™ are also 
considered in [T31 HH HBl HZl UHl 1201 [22] . 

In this paper, we propose a new approach to exponential estimates for solu- 
tions of systems of partial differential equations. Our approach is based on the 
limit operators method. This method was employed earlier to study the essential 
spectrum of perturbed pseudodifferential operators, which has found applica- 
tions to electro-magnetic Schrodinger operators, square-root Klein-Gordon, and 
Dirac operators under very general assumptions for the magnetic and electric 
potentials at infinity. Based on the limit operators method, a simple and trans- 
parent proof of the well known Hunziker, van Winter, Zjislin theorem (HWZ- 
Theorem) for multi-particle Hamiltonians was derived in [22, 23]. In [25], the 
limit operators method was applied to study the essential spectrum of discrete 
Schrodinger operators. 

Let 

A(x,D)u(x)= a a (x)D a u{x), x e R" 

\a\<m 

be a uniformly elliptic system of partial differential operators of order m on R™ 
with bounded and uniformly continuous N x N matrix- valued coefficients a a . 
Further let w '■= expv be a weight on R n such that lim^^oo w(x) = oo and 
assume that all first and higher order derivatives of v exist and are bounded. 

With the operator A(x, D) and the weight w, we associate a family of partial 
differential operators, 

B t := A(x,D + UVv(x)), te [-1,1]. 

We let sp ess B t denote the essential spectrum of B t considered as an unbounded 
closed operator on the Hilbert space L 2 (R n ,C N ) with domain H m (R n ,C N ). 
Further we let L 2 QH n ,C N ,w) and H m (W n ,C N ,w) refer to the corresponding 
weighted spaces of vector- valued functions u such that wu € L 2 (W l ,C N ) and 
wu 6 H m (R n , C N ), respectively. The following result, which will be proved in 
this paper, provides exponential estimates for the decay of the solutions of the 
equation A(x, D)u = f. 

Theorem 1 Let A(x, D) be a uniformly elliptic on R n matrix partial differential 
operator with bounded and uniformly continuous coefficients, and assume that 
i sp ess B t for every t e [-1,1]. Let u € H rn {R n ,C N ,1V- 1 ) be a solution of 
the equation A(x, D)u = f with f € L 2 (R n , C N , w) . Then u € H m (R n , C N ', w). 

By this result, the derivation of exponential estimates of solutions of the equa- 
tion A(x, D)u = f is basically reduced to the calculation of the essential spec- 
trum of the Bt . The latter can be done by means of the limit operators method 
(see [26] [27] and the monograph [28]). For, one associates with each operator Bt 
the family op (St) of its limit operators B g t which, roughly speaking, describe 
the behaviour of the operator at infinity. Then it follows from [25] [27] [2S] that 

sp ess -B*= IJ sp-Bf, (1) 
s t 9 eo P (B t ) 
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where spBf denotes the spectrum of the unbounded operator Bf acting on 
L 2 (R", C N ). 

In many important instances, the limit operators are of an enough simple 
structure, namely partial differential operators with constant coefficients. Then 
formula fT} provides an effective tool to calculate the essential spectra of par- 
tial differential operators and thus, in view of Theorem [JJ to get exponential 
estimates for the solutions of the equation A(x, D)u = f. We will illustrate this 
statement by applying Theorem [JJ to verify the exponential decay and to obtain 
explicit exponential estimates for the eigenvectors of Schrodinger operators with 
matrix potentials, of Moisil-Theodorescu quaternionic operators with variable 
coefficients, and of Dirac operators. 

The contents of the paper is as follows. In Section 2 we introduce the main 
definitions and prove Theorem [TJ The following sections are devoted to several 
applications of Theorem [TJ We start in Section 3 with the essential spectrum 
and exponential estimates for eigenvectors of Schrodinger operators with matrix 
potentials. Note that the well known Pauli operator is a Schrodinger operator 
of this kind (see, for instance, [3]). Moreover, such operators appear in the 
Born-Oppenheimer approximation for polyatomic molecules jTOj [T5l [19]. For 
potentials which are slowly oscillating at infinity we describe the location of the 
essential spectrum and give exact estimates of the behavior of eigenfunctions of 
the discrete spectrum at infinity. 

In Section 4 we consider general Moisil-Theodorescu (quaternionic) operators 
with variable coefficients. Note that numerous important systems of partial 
differential operators of quantum mechanics, elasticity theory, and field theory 
admit a formulation in terms of quaternionic operators (see [3 [9J [UJ [12l [29] 
and the references cited there). We shall verify explicit necessary and sufficient 
conditions for quaternionic operators with variable coefficients to be Fredholm 
operators and derive exponential estimates at infinity for solutions of Fredholm 
quaternionic equations. 

In the concluding Section 5 we consider the Dirac operator on K 3 , equipped 
with a Riemannian metric, with electric and magnetic potentials which are 
slowly oscillating at infinity. 

2 Exponential estimates of solutions of systems 
of partial differential equations 

2.1 Essential spectrum 

We will use the following standard notations. 

• Given Banach spaces X,Y, C(X,Y) is the space of all bounded linear 
operators from X into Y. We abbreviate C(X, X) to C(X). 

• L 2 (M. n ,C N ) is the Hilbert space of all measurable functions on M" with 



3 



values in C N , provided with the norm 

1/2 



IMIl2(r~,c«) : = (J \\u(x)\\%. N dx 



• The unitary operator V h of shift by h £ R™ acts on L 2 (R n , C N ) via 
(Vfcu)(x) := u(x - h). 

• Ci,(R™) is the C*-algebra of all bounded continuous functions on R". 

• C^(R n ) is the C*-subalgebra of Cb(R n ) of all uniformly continuous func- 
tions. 

• SO(R n ) is the C*-subalgebra of C^(R") which consists of all functions a 
which are slowly oscillating in the sense that 

lim sup \a(x + y) — a(x)\ = 

for every compact subset K of R" . 

S0 1 ( 
that 



• S0 1 (M. n ) is the set of all bounded differentiable functions a on R™ such 



lim = Ofor j = 1, 

x^oo OXj 



.1 



Evidently, S0 1 (M n ) C 50(R"). 

We also use the standard multi-index notation. Thus, a = (a\, a n ) with 
aj e N U {0} is a multi-index, \a\ = ai + . . . + ot n is its length, and 

d a :=d^...d^ and L> a := H^)" 1 . . . (-id Xn ) a " 

are the operators of a th derivative. Finally, {£) :— (1 + l^l 2 ) 1 ^ 2 for £ G R n . 
We consider matrix partial differential operators of order m of the form 

(Au)(x)= a a (x)(D a u){x), x <E R™ (2) 

\a\<m 

under the assumption that the coefficients a a belong to the space 

C 6 "(R", C(C N )) := C 6 u (R n ) ® £(C N ). 

The operator A in (J2]) is considered as a bounded linear operator from the 
Sobolev space iJ m (R™,C") to L 2 (R n ,C N ). The operator A is said to be uni- 
formly elliptic on R" if 



inf 



N 

det aQ,(x)oj Q 

\a\—m 



>0 (3) 
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where S™ -1 refers to the unit sphere in R™. 

The Fredholm properties of the operator A can be expressed in terms of its 
limit operators which are defined as follows. Let h : N — > R n be a sequence 
which tends to infinity. The Arcela-Ascoli theorem combined with a Cantor 
diagonal argument implies that there exists a subsequence g of h such that the 
sequences of the functions x i— > a a (x + g(k)) converges as k — > oo to a limit 
function a£ uniformly on every compact set K C K n for every multi-index a. 
The operator 

A* := J2 < D<X 

a|<m 

is called the limit operator of A defined by the sequence g. Equivalently, A 9 
is the limit operator of A with respect to g if and only if, for every function 

xeC °°(K"), 

in the space C(H m (W l , C N ), L 2 (W l , C N )) and 

lim V_ g{m) A*V g{m)X I = (A g )* X I 
in C(L 2 {M n 1 C N ),H-' m (R n ,C N )). Here, 

A*u = ^ D a (a* a u) and (A 9 )*u = ^ D a {(a 9 a )*u) 

\a\<m \a\<m 

refer to the adjoint operators of A, A g : iJ m (R", C^) -> L 2 (R", C N ). Let finally 
op (A) denote the set of all limit operators of A obtained in this way. 

Theorem 2 Let A be a uniformly elliptic differential operator of the form 
(0). Then A : H rn (« n ,C N ) -> L 2 (R n ,C N ) is a Fredholm operator if and 
only if all limit operators of A are invertible as operators from H m (R n ,C N ) 
to L 2 (R n ,C N ). 

Proof. Let (I — A) a be the pseudodifferential operator with symbol (1 + |£| 2 ) Q > 
and let A be a partial differential operator of the form @. Then the operator 
B := A(I — A)~ m / 2 belongs to the matrix Wiener algebra considered in [21]. It 
follows from results of [23] that B, considered as an operator from L 2 (R n ,C N ) 
to L 2 (M. n ,C N ), is a Fredholm operator if and only if all limit operators of B 
are invertible on L 2 (M. n ,C N ). Since the limit operators B g of B and A g of A 
are related by the relation B g = A g (l — A)~" 1 / 2 (which comes from the shift 
invariance of the operator A), the assertion follows. ■ 

The uniform ellipticity of the operator A implies the a priori estimate 

IM|j?'»(K»,c* r ) < C (II^w||l2(B",c«) + IMIl^b-m:")) ■ ( 4 ) 

This estimate allows one to consider the uniformly elliptic differential operator A 
as a closed unbounded operator on L 2 (R™, C N ) with dense domain H m (W l , C ). 
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It turns out (see [3]) that A, considered as an unbounded operator in this way, 
is a (unbounded) Fredholm operator if and only if A, considered as a bounded 
operator from H m (W\C N ) to L 2 (W\C N ), is a (common bounded) Fredholm 
operator. 

We say that A G C belongs to the essential spectrum of A if the operator 
A — XI is not Fredholm as an unbounded differential operator. As above, we 
denote the essential spectrum of A by sp ess A and the common spectrum of A 
(considered as as unbounded operator) by sp A. Then the assertion of Theorem 
[2] can be stated as follows. 

Theorem 3 Let A be a uniformly elliptic differential operator of the form (0) • 
Then 

s P ess ^= IJ SPA,- (5) 

A g Gop A 

2.2 Exponential estimates 

Let it; be a positive measurable function on R", which we call a weight. By 
L 2 (M. n , C N , w) we denote the space of all measurable functions on 1" such that 

IM|l2(r~ ,€«,«>) : = IMIl2(r™,c«) < oo. 

In what follows we consider weights of the form w = exp v where d Xj v € (R™ ) 
for j = 1 , . . . , n and 

lim d 2 x .v{x) = for 1 < i, j < n. (6) 

x — >oo 2 3 

We call weights with these properties slowly oscillating and let 1Z stand for the 
class of all slowly oscillating weights. 

Examples of slowly oscillating weights can be constructed as follows. Given 
a positive C°°-function I : 5 n_1 — * M, set Vi(x) := Z(x/|x|)|x|. Then u>i := expvi 
defines a weight on R™. Clearly, vi is a positively homogeneous function, that 
is vi(tx) = tvi(x) for alH > and x £ MP. Moreover, v t £ C°°(R™ \ {0}), and 
Wvi(oj) = l(u)u) for every point u> £ S 1 ™ -1 . Let -D; refer to a C°°-function on 1" 
which coincides with vi outside a small neighborhood of the origin. Then the 
weight u>i := expu; belongs to the class TZ. Moreover, 

lim Vvi(x) = Vvi(uj) = 1(u)uj (7) 

for lu £ S n -\ 

Proposition 4 Let A be a differential operator of the form (0), and let w — 
exp v be a weight in TZ. Then 

w~ 1 Aw= a a {D + iVv) a + <S> + R 

\a\<m 

where R := Yl\ a \<m-i baD a is a differential operator with continuous coeffi- 
cients such that lim-r^oo b a {x) = 0. 
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For a proof see [H] where a similar result is derived for a large class of pseudo- 
differential operators. The following is taken from [8], p. 308. 

Proposition 5 Let X\. X%, Y\ and Y<i be Banach spaces such that X\ is densely 
embedded into X2 and Y\ is embedded into Y%. Further let A : X2 — > I2 and 
A\x t '■ X\ — > Y\ be Fredholm operators, and suppose that 

ind (A: X 2 —> Y 2 ) = ind (A\ Xl : X l -> Y x ). 

If u £ X 2 is a solution of the equation Au — f with right-hand side f £ Y\, then 
u G Xi. 

Theorem 6 Let A be a uniformly elliptic differential operator of the form (0) , 
and let w — exp-u be a weight in 1Z such that lim^^oc w(x) = +00. For t G 
[—1, 1], set 

A w , t := a a (D + iWv) a , 

\a\<m 

and assume that 

0^ |J s Pe«A 0> t= (J U S P^' ( 8 ) 

te[-i,i] «e[-i,i]A» it eop(A w , t ) 

If u is a function in H m (R n , C , w^ 1 ) for which Au is in L 2 (M. n , C , w), then 
u already belongs to H m (R n , C , w). 

Proof. Note that A : H m (W\ C N , to*) -> L 2 (W\C N ,w t ) is a Fredholm opera- 
tor if and only if w^Aw 1 : H m (R n , C N ) -> L 2 (R n , C N ) is a Fredholm operator, 
and that the Fredholm indices of these operator coincide. Proposition [4] implies 
that 

w^Aw* = a a (D + itVv) a + R t , 

\a\<m 

where R t = J2 \ a \ <m-i b a ,tD a and lini x _ >O0 b a ,t(x) = for every t G [—1,1]. 
Hence, op(w~ t Aw t ) = o~p(A W}t )- Moreover, it is not hard to see that the 
coefficients b a t depend continuously on t G [—1, 1]. Hence, the family w^ t Aw t 
depends continuously on t, and condition (JSJ implies that all operators w~ t Aw t : 
H m (R n ,C N ) -> L 2 (R n ,C N ) with t G [-1,1] are Fredholm operators and that 
the Fredholm indices of these operators coincide. This implies that each operator 

A : H m (R n ,C N ,w t ) -> L 2 (R n ,C N ,w t ) 

owns the Fredholm property and that the index of this operator is independent of 
t G [-1, 1]. Since H m (R n , C N , w) is densely embedded into H m (R n , C N , to" 1 ), 
we can apply Proposition [5] to obtain that all solutions of the equation Au = f 
with right-hand side / G L 2 (R™, C", to), which a priori are in H rn (R n , C N , w^ 1 ), 
in fact belong to H m (R n , C N , w). m 
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Corollary 7 Let A be a uniformly elliptic differential operator of the form (0), 
and let w = expv be a weight in 1Z with lim^^oo w(x) — +00. Let A £ sj> dis A 
and A ^ sp ess A tw for all t £ [0, 1] . Then every eigenfunction of A associated 
with A belongs to the space H m (M. n , C N , w). 

Indeed, this is an immediate consequence of Theorem [6] since eigenfunctions 
of uniformly elliptic operator of order m necessarily belong to H m (M. n , C N ). 

3 Schrodinger operators with matrix potentials 

3.1 Essential spectrum 

We consider the Schrodinger operator 

H := (id Xj - a 3 )p> k {id Xk - a k )E + $ (9) 

where E is the N x N unit matrix, a — (a%, ...,a„) is referred to as the the 
magnetic potential, and $ = ($ pq )p q= i is a matrix potential on R™, the latter 
equipped with a Riemann metric p = (pjk)™h=i wmcn is subject to the positivity 
condition 

inf p ik (x)<Jw k > 0, (10) 
lei-sues™- 1 

where Pjk(x) refers to the matrix inverse to pP k (x). Here and in what follows, 
we make use of Einstein's summation convention. 

In what follows we suppose that fP k and aj are real-valued functions in 
S , 1 (M") and that Q pq £ SO(M. n ). Under these conditions, H can be considered 
as a closed unbounded operator on L 2 (R n , C N ) with domain H 2 (R n , C N ). If $ 
is a Hermitian matrix-valued function, then Tl is a self-adjoint operator. 

The limit operators of Ji are the operators with constant coefficients 

W = (id Xj -a3)4 k (id Xk -al)E + & 

where 

a 9 := lim a(g m ), p g := lim p(g m ), $ 5 := lira $(ff m ). (11) 

m — ^00 m — *oo m — ^oo 

The operator H. 9 is unitarily equivalent to the operator 

H{ :=-pi k d Xj d Xk E + ^, 

which on its hand is unitarily equivalent to the operator Hf of multiplication 
by the matrix-function 



8 



acting on L 1 (R™ , C ) . Evidently, 

s P n{ = (J r 
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where T 9 := fj, 9 + K and the |Lt|, 1 < j < N, run through the eigenvalues of the 
matrix $ 9 . Thus, specifying (J5J) to the present context we obtain the following. 

Theorem 8 The essential spectrum of the Schrddinger operator TL is given by 

N 
9 3 = 1 

where the union is taken with respect to all sequences g for which the limits in 
(f77j) exist. 

The description (fH?)) of the essential spectrum becomes much simpler if $ is a 
Hermitian matrix function, in which case TL is a self-adjoint operator. 

Theorem 9 Let the potential <£> be a Hermitian and slowly oscillating matrix 
function. Then 

sPess ft = M*,+oo) 

where 

<i$ := liminf inf (<&(x)(p,tp). 

x— >oo ||<p|| = l 

Proof. Since <£> 9 is Hermitian matrix, 

7 9 := inf ($ g ip, ip) 
ll¥>ll=l 

is the smallest eigenvalue of $ 9 . Hence, spTL 9 — [7 ff ,+oo) and, according to 

SPess K = (J ^9' +°°) = t illf ^9' +°°)- 
9 

It remains to show that 

inf 7„ = d$. (13) 

9 

Let g be a sequence tending to infinity for which the limit 

:= lim $(g(m)) 

m— ' — 

exists. Then, for each unit vector tp £ 

($ 9 tp,(p) = lim ($(g(m))ip,(p) > liminf (&(g(m))<p, ip) > d$, 

m — >oo x — >oo 

whence j g > g?$. For the reverse inequality, note that there exist a sequence 
go tending to infinity and a sequence ip in the unit sphere in C N with limit p>o 
such that 

d$ = lim (®(go(m))(p m ,y m ) = (^Vo^o) > 7 S o- 

m — >oo 

Thus, 7 9o = o?$, whence (fTB")) . ■ 



m — >oo 
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3.2 Exponential estimates of eigenfunctions of the discrete 
spectrum 

Here we suppose that the components p 1 of the Riemann metric, the coefficients 
a a and the weight w are slowly oscillating functions and that $ is a Hermitian 
slowly oscillating matrix function. Every limit operator (w~ 1 H.w) g of w^Tiw 
is unitarily equivalent to the operator 

n g w := ff(D Xj + i(Vv)<>)(D Xk + i(Vv) 9 k ) + 

where p^ g k and $ s are the limits defined by (fTTj) and 

(Vv) 9 := lim (Vv)(g(k)) G K™. (14) 

— ► OO 

We set 

|(V V )|2:=^(V^(V«) fe and |N^ s :=^(Vt,)f(V,)f. 

The operator 7ig, is unitarily equivalent to the operator of multiplication by the 
matrix-valued function 

W?(0 := + i(yv)l )(& + i(Vv) 9 k ) + i G M", 

the real part of which is 

) = tftek + ($ s - \(Vv)% g E). (15) 
Corollary implies the following. 

Theorem 10 Let A G sp dis 7i, and let w = expv be a weight in 1Z for which 



lim sup | (Vv) (x) | p ( x -) < y/d$ - A. 

Then every X-eigenf unction ofTi belongs to H 2 (M n 1 w). 
Corollary 11 Let A G sp dis H, and let c € M satisfy 



< c < 

psup 

where 

p BU P := liminf sup (p»>(x)u j Uk) 1/2 ' 

XTien i/ie every \-eigenf unction ofH belongs to the space H 2 (M. n ,C N ,w) with 
weight w(x) — e c ^ . 
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4 Quaternionic operators 



We let H(C) denote the complex quaternionic algebra, which is the associa- 
tive algebra over the field C of complex numbers generated by four elements 
1, ei, e 2 , e 3 subject to the conditions 

eie 2 = e 3 , e 2 e 3 = ei, e 3 ei = e 2 

and 

l 2 = 1, e 2 , = -1, le k = e k l = e k , e 3 e k = -e k ej 

for j, k = 1,2, 3. Each of the elements 1, e\, e 2 , e 3 commutes with the imaginary 
unit i. Hence, every element q € H(C) has a unique decomposition 

q = 9o + qiei + q 2 e 2 + g 3 e 3 =: q Q + q 

with complex numbers qj . The number go is called the scalar part of the quater- 
nion q, and q is its vector part. One can also think of H(C) as a complex linear 
space of dimension four with usual linear operations. With respect to the base 
{1, ei, e 2 , e 3 } of this space, the operator of multiplication by 1 has the unit ma- 
trix E4 as its matrix representation, whereas the matrix representations jj of 
the operators of multiplication by ej, j — 1,2,3, are real and skew-symmetric, 
that is 7] = — 7j. The space H(C) carries also the structure of a complex Hilbert 
space via the scalar product 

(q, r}w(c) ■= qon) + qrn + 92^2 + q^- 

By L 2 (]R 3 , H(C)) we denote the Hilbert space of all measurable and squared 
integrable quaternion- valued functions on M 3 which is provided with the scalar 
product 

(u,v) L 2 {R 3 Mcy) := / (u(x),v(x))m C )dx. 

In a similar way, we introduce the quaternionic Sobolev space 7? 1 (R 3 , H(C)). 
Further we write M v for the operator of multiplication from the right by the 
complex quaternionic function ip, that is 

(M v u)(x) = u(x)<p(x) for x e M 3 . 

Clearly, if ip G L°°(R 3 , H(C)), then M v acts as a bounded linear operator on 
L 2 (M 3 ,H(C)). 

Differential operators of the form 

3 

A(x,D)u(x) ■.= ^2aj(x)D Xj e j u(x) + M' p( -^u(x), x S M 3 , (16) 
3=1 

can be considered as generalized Moisil-Theodorescu operators. Note that each 
operator of the form (fT6|) corresponds to a matrix operator with respect to the 
basis {1, ei, e 2 , e 3 }. It has been pointed out in [HI HU [H] that some of the 



11 



most popular operators of mathematical physics, including Dirac and Maxwell 
operators, are of the form (fl~6)) . 

In this section, we suppose that the coefficients a,j belong to S'0 1 (R 3 ) and 
satisfy 



inf \aj(x)\ > for j = 1,2,3 



(17) 



and that the components ipk of ip belong to .SC^R 3 ) 
The main symbol of the operator A is 



A)(x,0 = ^a j {x){iS, j )e J 
3=1 



Hence, 



i=i 

is a scalar function, and from (| 1 T[) we conclude that the associated operator Aq 
is uniformly elliptic on R 3 . 

Theorem 12 The quaternionic operator A(x, D) thought of as acting from 
7J 1 (R 3 ,H(C)) to L 2 (R 3 ,H(C)) is a Fredholm operator if and only if 



lim inf 

x — >oo 



4>(s,0+£vj(aO S 
3=1 



> for every £ G 



(18) 



Proof. The limit operators of A(x, D) are the operators with constant coeffi- 
cients 

3 

A 9 (D)u := ojA ' + M<p9u - 

3=1 

Let A»(£>) := tfD^.ej - M* 9 . Then 

A^p)i«(£>) = -^(af) 2 ^.-(^) 2 

3=1 

where 

-(^) a = (^) 2 + (^) 2 + (^) 2 . 

Condition IJI7JI implies that A9(D) : ii^R 3 , H(C)) -> L 2 (R 3 , H(C)) is an in- 
vertible operator if and only if 



inf 



3=1 



> for every £ e 



(19) 



Hence, all limit operators A 9 (D) of A(D) are invertible as operators from 
7J 1 (R 3 ,H(C)) to L 2 (R 3 ,H(C)) if and only if condition $TE$ holds. ■ 
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Theorem 13 In addition to the above assumptions, let all junctions <Zj and tpj 
be real-valued, and let w = expi> be a weight in 1Z with lim^^oo v(x) = +00. If 
the condition 

Hmmf^^ W -«l(-)(^) 2 j >0 (20) 

is satisfied, then every solution u £ i? 1 (R 3 ,H(C), of the equation Au = f 

with right-hand side f G L 2 (R 3 , H(C), w) belongs to H 1 (R 3 , H(C), w). 

Proof. Let t £ [—1,1]. The limit operators of A w t(x,D) are operators with 
constant coefficients of the form 

Al, t (D) =Y / a a j (D Xj +it(^-y)e j + MV 9 . 



j 



As above, let 



dv 



A 9 w , t (D) :=Y f a°(D„ j +it{—)°)e j -M*'. 

j=l 3 

Then A 9 W t (D)A 9 w t (D) is a scalar operator with symbol 

KASHAS = EK) 2 te + i^) 9 f + E^-) 2 ^ 

the real part of which is 

3 3 

^(,n?\ 2 -t 2 \l 



w« it (0<t(0) - E«) 2 £f + E(^f) 2 - * 2 [(^ L ) S ] 2 - 



i=i j'=i 
Condition (|20f implies that 

^,t(0<*(0^0 (21) 

for every £ G R 3 and t £ [—1,1]. Without change of notation, we now consider 
^«;t(0 as a 4 x 4 matrix- valued function. The matrix A 9 W 4 (£) is invertible 
for every £ G M 3 and t £ [-1,1] and for every sequence g which defines a 
limit operator. Together with condition (fT7|) . this fact implies that A 9 W t (D) : 
H^R 3 ,!!^)) -> L 2 (R 3 ,H(C)) is an invertible operator for every t £ '[-1,1] 
and for every sequence g which defines a limit operator. Hence, Theorem 1131 is 
a consequence of Corollary ■ 



5 Dirac operators 

5.1 Essential spectrum of Dirac operators 

In this section we consider the Dirac operator on R 3 equipped with the Riemann 
metric tensor (pjk) depending on x G R 3 (see for instance [30 ) . We suppose 
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that there is a constant C > such that 

Pjk(x)€ j € k > C|C| 2 for every x G R 3 (22) 

where we use the Einstein summation convention again. Let p* k be the tensor 
inverse to pj k , and let (j>> k be the positive square root of p> k . The Dirac operator 
on R 3 is the operator 

V ■= pk(^ k Pj + Pj^ k ) + loc 2 m - e$ (23) 

acting on functions on R 3 with values in C 4 . In (|23[> . the j k , k = 0, 1, 2, 3, are 
the 4x4 Dirac matrices, i.e., they satisfy 

Jjjk + Ikl] = 26 jk E (24) 

for all choices of j, k = 0, 1, 2, 3 where E 1 is the 4x4 unit matrix, 

». ' 77^7 J = 1.2, 3, 

where /i is the Planck constant, A = (Ax, A2, A3) is the vector potential of the 
magnetic field H, that is H = rot A, $ is the scalar potential of the electric field 
E, that is E = grad $, and m and e are the mass and the charge of the electron. 

We suppose that p> k , Aj and 4> are real- valued functions which satisfy the 
conditions 

p> k e S'O^R 3 ), ^ e SO^R 3 ), $e50(M 3 ) (25) 

for j, k — 1,2, 3. We consider the operator T> as an unbounded operator on the 
Hilbert space L 2 (R 3 ,C 4 ) with domain iJ^R^C 4 ). The conditions imposed on 
the magnetic and electric potentials A and $ guarantee the self-adjointness of 
T>. The main symbol of T> is 

ao(x,0 = c<j) jk (x)^jk- 

Using (|24| and the identity (j>> k (x)(j) rt (x)5 k t — p' r ( x ) we obtain 

a (x,Z) 2 = c 2 h 2 ^ k {xW\x)^ rlklt 
= c 2 h 2 4P k {x)4> rt {x)5 kt ^ r 
= c 2 h 2 p> r (x)^C r E. (26) 

Together with (|22[) this equality shows that T> is a uniformly elliptic differential 
operator on R 3 . 

Conditions (|25|) imply that limit operators T> g of T> defined by the sequences 
g : N — > Z 3 tending to infinity are operators with constant coefficients of the 
form 

^ 9 - c 7fe 0f (Dj + -Af) + 7o mc 2 - (27) 
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where 

4P a k := lim 4P k {g{m)), A 3 , := lim AJgtm)), <I> 9 := lim $fa(m)). (28) 

Note that one can suppose without loss of generality that the sequence g is 
such that the limits in (f2"8")) exist. In the opposite case we pass to a suitable 
subsequence of g. 

The operator T> g is unitarily equivalent to the operator 

V\ = cjtufDj + l0 mc 2 - e$ 9 , 

and the equivalence is realized by the unitary operator 

T Ag :f^e^ Aa -*f where A> := (M,^A%). 

Let 

$ sup := lim sup $ (x), $ inf := liminf 

a; — >oc ^ *°° 

Then the interval [$ , <£> sup ] is just the set of all partial limits $ s of function 
$ as a; — * oo. 

Theorem 14 Let conditions \25\ be fulfilled. Then the Dime operator 

V : H\R 3 ,C 4 ) -> L 2 (M. 3 ,C 4 ) 

is a Fredholm operator if and only if 

[$ inf , $ sup ] C (-mc 2 /e, mc 2 /e). (29) 

Proof. Set Dg(f) := ch lk ^ g k ^ + j mc 2 and := lim^^ ^' fe ( 3m ). Then 

(Hg(0-e$?S) (2?g(0 + e*»£?) 

= (c 2 h 2 pf^ k + m 2 c 4 - (e&) 2 ) E. (30) 

Let condition (|2"9"|) be fulfilled. Then every partial limit <f> 9 = lim^oo <&{g(k)) 
of $ lies in the interval (— mc 2 /e, mc 2 /e). The identity (|3H)l implies that 

det(P»(0-e$ 9 L;) t^O 

for every £ G tt 3 . Hence, the operator V] : iJ^R^C 4 ) -> L 2 (R 3 ,C 4 ) is invert- 
ible and, consequently, so is T> g . By Theorem [2l I? is a Fredholm operator. 

For the reverse implication, assume that condition (f2"9"| is not fulfilled. Then 
there exist a number <£> 9 ^ (— mc 2 /e, mc 2 /e) and a vector £° £ K 3 such that 

c 2 ^°4°+m 2 c 4 -(e^) 2 =0. 

Given £° we find a vector u e C 4 such that v := (X>g(^°) + (e^ 9 )^) u ^ 0. 
Then fl30|) implies that 

(23g(£°) - e$*£) « = 0, 

whence 

det(Pg(f°)-e$ 9 £;) = 0. 

Thus, the operator 2? g is not invertible. By Theorem[2l I? cannot be a Fredholm 
operator. ■ 
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Theorem 15 If condition t\25§ is satisfied, then 

s Vess V = (-00, -e$ inf - mc 2 ] U [-e$ sup + mc 2 , +00). 

Proof. Let A e C. The symbol of the operator T> g — XI is the function £ 1— ► 
^0(0 _ ( e$9 + X ) E - Invoking (001) we obtain 

(X>g(0 - (e$» + A)£) (2?g(0 + (e$» + A)£) 

= (cVf + m 2 c i - (e$ 9 + A) 2 ) E. (31) 

Repeating the arguments from the proof of Theorem 1 14| wc find that the eigen- 
values Aj_(£) of the matrix 2?q(£) — e<fr\E are given by 

A|(£) ■= -e$ 9 ± (c 2 4%Zk + m 2 c 4 )V2. (32) 

From |32]) we further conclude 

{A € R : A = Ai(£), £ e K 3 } = (-00, -e$ 9 - mc 2 ], 
{A € K : A = \ q + (£), £ E R 3 } = [-e$ 9 + mc 2 , +00). 

Hence, 

spV 9 = (-co, -e$ s - mc 2 ] U [-e$ s + mc 2 , +00), 
whence the assertion via Theorem [3J ■ 

Thus, if $ sup - $ inf > 2mc 2 /e, then sp ess V is all of R, whereas sp ess V has a 
proper gap in the opposite case. 

5.2 Exponential estimates of eigenfunctions of the Dirac 
operator 

Theorem 16 Let the conditions <\25\ be fulfilled. Let A be an eigenvalue ofT> 
which lies in the gap (— e$ mt — mc 2 , — e$ sup + mc 2 ) of the essential spectrum. 
Further, let w — expw be a weight in 1Z with linxr^oo w(x) = 00 which satisfies 



limsup \S/v(x)\ p[x) < —^/m?c^ - (e$ su P + A) 2 . (33) 

x — >oo Cll 

Then every eigenfunction ofD associated with A belongs to H 1 (M. 3 , C 4 , w). 

Proof. Let A e (-e$ mf -mc 2 , — e$ sup + mc 2 ) be an eigenvalue of T>. As above, 
we examine the spectra of the limit operators {T> w ,t)g of T> Wt t := w~ t T>w t for 
t running through [0, 1]. Let (T> w ,t)g be a limit operator of D w ,t with respect 
to a sequence g tending to infinity. One easily checks that (T> tw ) g is unitarily 
equivalent to the operator 

(V tw ) g := A t , g - e&E 
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where 

A t , g := c 7fc 0f (Dj + ith (-Q^y) + lomc 2 . 
The operator A t . g has constant coefficients, and its symbol is 



Further, 

m 



m 



ox k 



m [(m 2 c 4 - (e$ 9 + A) 2 ) E] 

dv dv 



(e$ 9 + A) 5 



E 



Assume that condition 



is fulfilled. Then, since c 2 h 2 p> g k £,j£,k > 



fcf.7„,«,A)>0 



for all t € [0, 1] and for all sequences g — > oo for which the limit operators 
exist. Hence, (|33p implies that the matrix A t — (e$ 9 — X)E is invertible for 
every (el 3 . On the other hand, due to the uniform ellipticity of A t>g one has 

A S sp \ A ti g(£) - (e$ 9 - A).Ej if and only if there exists a £ S M 3 such that 

the matrix A tjg ((;o) — (e& 9 + \)E is not invertible. Thus, A ^ sj)(V tw ) g for every 
t 6 [0, 1] and every sequence g — > oo. Via Corollary [3 the assertion follows. ■ 



We conclude by an example. Let the conditions (|25| be fulfilled, and let A 

-e$ sup + mc 2 ) and u\ an associated 



inf 



mc 



be an eigenvalue of T> in (— e<!> 
eigenfunction. If a satisfies the estimates 



< a < 



y /m 2 c i - (e$ su P + A) 2 



where 



chp su P 

9 sup := liminf sup ( i p> k (x)u j w k ) 1/2 , 



luES 2 



then u x S i7 1 (M 3 ,C 4 ,e a < ;r >). 
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